The Journal of Physical Chemistry Letters
● American Chemical Society (ACS)
Preprints posted in the last 90 days, ranked by how well they match The Journal of Physical Chemistry Letters's content profile, based on 63 papers previously published here. The average preprint has a 0.03% match score for this journal, so anything above that is already an above-average fit.
Hungerland, J.; Timmer, D.; Frederiksen, A.; Lünemann, D. C.; Thöle, D.; Saberamoli, G.; Schmidt, J.; Kumar, K.; Bartölke, R.; de Sio, A.; Mouritsen, H.; Lienau, C.; Solov'yov, I. A.
Show abstract
Photoactivated intermolecular electron transfer (ET) in cryptochromes proceeds along chains of aromatic residues and creates a spatially separated pair of radical electrons. Ultrafast time-dependent spectroscopy can provide experimental insight into this process and theoretical estimates of charge transfer rates are commonly obtained via Marcus theory. Here, we present a new perspective on the ET in European robin cryptochrome 4a (ErCry4a) that synthesizes insights from real-time ET calculations, ultrafast spectroscopic measurements and analytical derivations. The simulations exemplify that molecular vibrations play an essential role in enabling the ET dynamics, which was further rationalized through analytical derivations. Ultrafast pump-probe spectroscopy provided experimental access to the first 1.5 ns of the ET cascade, where multiple radical pair recombination rates arise due to the dynamic equilibrium along the ET chain. We show that the motions of the protein environment and the ET dynamics are inseparably coupled, violating the timescale separation required for Marcus theory. The presented results highlight that non-equilibrium coupling between electronic and nuclear motion dominates ET kinetics in ErCry4a during the first nanosecond after photo-excitation. The findings exemplify the limits of Marcus theory and refine the interpretation of ultrafast spectroscopic signatures in cryptochromes.
Kedjar, Y.; Hognon, C.; Douki, T.; Dumont, E.; MONARI, A.
Show abstract
The repair of photo-induced DNA lesions through nucleotide excision repair machinery is still the source of important questions. It has been observed that the repair rate of the different cyclobutane pyrimidine dimers, i.e. the photoproducts induced by dimerization of two {pi}-stacked pyrimidines (T<>T, T<>C, C<>T, C<>C), depends on the nucleobases involved in the lesion. TT derivatives (T<>T) are removed more slowly than those containing cytosine, especially in 5. Using all-atom molecular dynamics simulations and free-energy calculations, we demonstrate that the variation of the repair rate observed in human skin and in cultured cutaneous cell is associated to the recognition of the four lesions by the DDB2 protein moiety, and more specifically by the differential structural deformation induced on the complementary strand. Indeed, while C<>C and C<>T induce a larger deviation on the groove parameters, T<>T and T<>C, instead, affect DNA structure to a lesser extent. less affected. These effects then hamper differentially the downstream recruitment of the repair complexes. The observed DNA deformation correlates with the experimental repair rate and provides a structural rationale for the different repair rates of CPD by nucleotide excision repair machinery. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=105 SRC="FIGDIR/small/724087v1_ufig1.gif" ALT="Figure 1"> View larger version (43K): org.highwire.dtl.DTLVardef@cf6b6dorg.highwire.dtl.DTLVardef@195e35forg.highwire.dtl.DTLVardef@1829296org.highwire.dtl.DTLVardef@165baba_HPS_FORMAT_FIGEXP M_FIG C_FIG
Smyth, S.; Liu, Z. H.; Tsangaris, T.; Head-Gordon, T.; Forman-Kay, J. D.; Gradinaru, C. C.
Show abstract
Eukaryotic cap-dependent translation initiation is regulated by binding of the predominantly folded eukaryotic initiation factor 4E (eIF4E) to the intrinsically disordered eIF4E binding proteins (4E-BPs). Here, we report full-length atomistic conformational ensembles generated by IDPConformerGenerator and optimized by X-EISDv2 workflow for both apo 4E-BP2, the neuronal 4E-BP, and 4E-BP2 in complex with eIF4E, using data from single-molecule fluorescence and nuclear magnetic resonance (NMR), together with select coordinates from a 4E-BP1:eIF4E crystal structure. Structural sampling within dynamic complexes is often under-appreciated, with NMR and crystal structure data for 4E-BP:eIF4E suggesting different degrees of structural heterogeneity. Our ensemble models validated by solution spectroscopy data enable comparison of free 4E-BP2 and its complex with eIF4E. This shows a delocalization of contacts around canonical regions, which supports previous findings of unidirectional conditional occupancy of the binding sites. Two new contact regions emerged: one between the disordered N-termini of eIF4E and 4E-BP2, which may play an allosteric role in tuning the binding affinity, and the other between the C-terminus of 4E-BP2 and an extended region of eIF4E, which is consistent with the extended, dynamic binding interface that we reported previously. These results support a model of translation regulation in which the dynamic 4E-BP2:eIF4E complex facilitates accessibility of regulatory sites of 4E-BP2 when bound.
Wieners, L.; Garcia, M. E.
Show abstract
Ultraviolet (UV) radiation induces DNA damage associated with cancer and aging, yet the sequence dependence of UV absorption remains to be investigated. Here, we present a systematic study of the UV absorption spectra of DNA based on all-electron Hartree-Fock calculations. We analyze all possible sequences up to four base pairs, as well as longer randomized sequences and genomic nullomers - motifs which are missing in a given genome. We observe a pronounced sequence dependence: cytosine- and guanine-rich motifs exhibit significantly enhanced absorption, whereas adenine-thymine-rich sequences absorb up to four times less in the mid-UV range. Notably, the human genome is biased toward adenine-thymine-rich sequences, giving it an increased susceptibility to UV-induced damage. In addition, we introduce a computational framework enabling spectral calculations of large DNA and RNA fragments, opening the door to large-scale optical analyses.
Marciniak, A.; Kozielewicz, P.; Mitrovic, D.; Delemotte, L.
Show abstract
Cells communicate with their environment by integrating signals, often chemical in nature, triggered by specific molecules bind to specific membrane-bound receptors, resulting in a downstream signaling cascade. Arguably, G-protein-coupled receptors (GPCRs) constitute the most pharmacologically important family of such receptors, binding small molecules, peptides, lipids, and hormones with high specificity. However, despite a highly conserved fold and sequence similarity, GPCRs are still mostly studied on a case-by-case basis. Here, we infer a general, evolutionarily conserved mechanism of class A GPCR activation. By leveraging coevolution and machine learning methods applied to all class A GPCRs structures, we derive a mathematical description (a so-called collective variable - CV) of the receptor's activation state which is independent of its sequence. Then, we bias molecular dynamics simulations along this CV to obtain transitions between activation states of a diverse set of class A GPCR family members. To demonstrate that our model generalizes beyond GPCRs in our training set, we obtain conformational transitions of an orphan receptor, GPR183. Finally, we show that we can model ligand effect on the receptors by converging Free Energy Surfaces of activation of the {beta}2-adrenergic receptor within this common mechanism framework. These results, to our knowledge, prove for the first time the existence of a mechanism uniting all class A GPCRs. Our approach thus facilitates direct comparisons between receptors and opens up the possibility of structural and dynamical studies of many orphan and understudied GPCRs. It also serves as a blueprint for inferring family-wide protein mechanisms.
Tsugawa, S.; Kikuchi, K.; Date, K.; Nonoyama, T.; Kang, Z.; Ueno, T.
Show abstract
Spiral geometries commonly occur in natural and engineered systems and are fundamentally described by curvature and torsion. In deformation-dominated systems, these variables evolve dynamically, requiring a continuum mechanical framework to link geometry and deformation. This study focused on refractile bodies (R-bodies), protein supramolecular assemblies that undergo reversible roll-spiral transformations in response to stimuli such as pH changes. Although multiple R-body types with distinct morphologies and unrolling behaviours were experimentally identified, their deformation mechanisms lack quantitative theoretical descriptions. We proposed a deformation-gradient-tensor-based continuum model incorporating geometrical mapping from the rolled to spiral state within a unified framework. The model successfully reconstructed macroscopic deformation behaviours of types 51, 7, and Pa R-bodies by capturing differences in unrolling behaviours, tapered geometry, and spatio-temporal evolution. The analysis revealed that deformation proceeds through a coupled process in which the curvature decreases via straightening, while the torsion increases by twisting. Importantly, the framework connected the macroscopic morphology with microscopic lattice deformation, enabling quantitative inference of lattice intervals and angles. The proposed comprehensive geometric model of the R-body roll-spiral transformation offers a general mathematical foundation for understanding deformation-driven spiral transformations in soft matter systems.
von Roten, V.; Ivanovic, M. T.; Gopi, S. R.; Holla, A.; Prestel, A.; Nüesch, M.; Tamburrini, K. C.; Nettels, D.; Kragelund, B. B.; Best, R.; Schuler, B.
Show abstract
The free-energy surfaces that underlie the conformational distributions of intrinsically disordered proteins (IDPs) are shallow and lack the deep minima characteristic of stable, folded structures. However, even in the absence of secondary or tertiary structure, sequence patterning can lead to conformational preferences and changes in chain dimensions as a function of solution conditions. While patterning effects have received extensive attention from simulation and theory, there is little corresponding data from experiment. Here we investigate the impact of charge patterning on chain dimensions and dynamics in a set of specifically designed polyampholytic IDP variants across the natural range of charge segregation with single-molecule FRET, nanosecond fluorescence correlation, circular dichroism, and NMR spectroscopy. We find that the conformational ensembles and their cooperative response to salt concentration show prominent and systematic dependencies on charge patterning, and to some extent on residue type. In contrast, the chain dynamics remain in the tens-of-nanosecond range, consistent with the absence of pronounced free-energy barriers. In close combination with molecular simulations, we show how the concept of susceptibility can be used to quantify cooperativity in the absence of barriers and relate it to the shallow free-energy surfaces of IDPs.
Samajdar, R.; Chhabra, H.; Meigooni, M.; Yi, S.; Liu, X.; Wu, J. L.; Tajkhorshid, E.; Schroeder, C. M.
Show abstract
Stereochemistry underlies structure-function relationships across biology and materials science, ranging from proteins to electronic and spintronic materials. In this work, we investigate the electron transport properties of different oligopeptide stereoisomers using experiments and computational modeling. Single-molecule electronic experiments show that stereochemical modifications in tyrosine-based peptides lead to significant variations in molecular conductance along the peptide backbone due to enhanced stacking interactions and electronic coupling of aromatic side chains. In addition, stereochemical variations in alanine-based peptides give rise to changes in conductivity due to secondary structure interactions arising from {beta}-turn conformations. All-atom molecular dynamics (MD) simulations and quantum mechanical calculations are used to understand the molecular origins of the effect of stereochemistry on the structural and electronic properties of peptides. Overall, this work shows that stereochemical modification of non-terminal amino acids effectively controls electron transport due to aromatic side chain interactions or secondary structure effects. These insights open new avenues for the molecular design of peptide-based electronic materials with enhanced function.
Pradhan, S.; Tripathi, S. M.; Sharma, S.; Singh, A. P.; Sundriyal, S.; Patra, S.
Show abstract
G-quadruplex (GQ) structures within the HIV-1 long terminal repeat (LTR) regulate viral transcription and represent promising antiviral targets; however, detailed mechanistic understanding of their ligand recognition at the molecular level remains limited and has largely been investigated under dilute conditions despite the crowded and compartmentalized nature of intracellular environment. Here, we investigate the interaction of the cationic porphyrin TMPyP4 with the HIV-1 LTR-III GQ under dilute conditions and inside protein-rich phase-separated condensates that mimic intracellular biocondensates. Steady-state and time-resolved fluorescence measurements reveal a dual binding behavior that is not discernible from absorption spectroscopy. A high-affinity guanine-rich binding mode leads to efficient fluorescence quenching through electron transfer from ground-state guanine to excited TMPyP4, whereas a weaker non-guanine binding mode gives rise to enhanced and long-lived emission. Nucleotide-specific control experiments validate the origin of these distinct binding environments. Molecular docking and molecular dynamics simulations further support preferential binding of TMPyP4 at the terminal G-quartet together with a secondary binding mode near the quadruplex-duplex junction. Importantly, both TMPyP4 and LTR-III GQ preferentially partition into the condensates, where the hybrid GQ structure, dual binding behavior, and associated excited-state signatures remain preserved despite the crowded and viscous environment. Although a slight reduction in binding affinity is observed inside the condensates, the overall binding mechanism remains largely preserved due to compensatory effects arising from the condensate microenvironment. Overall, this work demonstrates that ligand recognition of viral GQ remains preserved within protein condensates and establishes fluorescence spectroscopy as a sensitive approach for resolving hidden binding heterogeneity in GQ-ligand interactions.
Kimura, T.; Katayama, T.; Ishikawa, S.; Mitra, S.; Yamano, Y.; Onizuka, K.; Nagatsugi, F.; Laha, S.; Naganathan, A. N.; Itoh, Y.; Takahashi, S.
Show abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) packages its single-stranded genomic RNA (gRNA) having about 30,000 nucleotides into virions by forming 35-40 granular ribonucleoprotein (RNP) units. Each RNP unit has a diameter of [~]15 nm. While it is generally assumed that the assembly of these RNPs is driven by the binding of the nucleocapsid (N) protein to the gRNA in the cytoplasm, the precise molecular mechanism remains to be fully elucidated. In this study, we develop an experimental strategy based on single-molecule fluorescence and fluorescence correlation spectroscopies to examine the formation of long-range base pairing within a candidate structural domain corresponding to nt 12230-12686 of the gRNA (gRNA12k). Our results demonstrate that the 5 and 3 regions of gRNA12k autonomously form long-range base pairing in near-physiological buffers containing mono- and divalent cations, independently of the N protein. This domain possesses an extensive secondary structure, is compact, and can unfold and refold reversibly upon heat treatment and cooling. Notably, the addition of the N protein melts the long-range base pairs, and causes the aggregation of multiple molecules of gRNA12k. Based on these observations, we propose a refined mechanism for the genome assembly in SARS-CoV-2: gRNA initially forms autonomous granular structures, which are subsequently reorganized and condensed by the N protein to chaperone the assembly of the entire gRNA. SignificanceSARS-CoV-2 organizes its exceptionally long genomic RNA (gRNA) having about 30,000 nt into 35-40 granular ribonucleoprotein (RNP) units for viral packaging. It has been assumed that the nucleocapsid (N) protein drives the formation of the RNP granules. In this study, we challenge this prevailing view by demonstrating that a specific region of the gRNA sequence inherently encodes the information to fold into a compact, granular architecture independently of any proteins. Unexpectedly, we found that the N protein partially melts the autonomous structures, suggesting that it acts as an RNA chaperone to facilitate flexible genome assembly. Our findings redefine the interplay between viral proteins and gRNA, offering a new perspective on the mechanism of coronavirus replication.
Feito, A.; Tejedor, A. R.; Ocana, A.; Teran, A.; Merlino, A.; Marasco, D.; Herrero, S.; R. Espinosa, J.
Show abstract
The inhibition of A{beta}42 ({beta}-amyloid) fibril formation is a key therapeutic strategy in Alzheimer's disease research. Paddlewheel diruthenium complexes have shown promising activity against A{beta}42 aggregation and preformed fibril disaggregation, yet their molecular mode of action remains poorly understood. In this work, we perform atomistic simulations to explore how charge modulation influences the interactions of three analogous paddlewheel diruthenium complexes, the parent neutral complex [Ru2Cl(D-p-FPhF)(O2CCH3)3], and its anionic [Ru2Cl2(D-p-FPhF)(O2CCH3)3]- and cationic [Ru2(D-p-FPhF)(O2CCH3)3]+ counterparts (D-p-FPhF- is the N,N' -bis(4-fluorophenyl)formamidinato ligand) with A{beta}42. Our results indicate that electrostatic tuning governs binding affinity and the extent of interaction across the A{beta}42 fibril surface. As the complexes' charge changes from -1 to +1, the interaction pattern shifts from localized contacts to widespread, multi-site engagement encompassing key charged, aromatic, and hydrophobic regions of A{beta}42. This enhanced binding correlates with longer-lived, thermodynamically stable interactions at the fibril interface, which effectively lower the free energy penalty for fibril disassembly. Overall, our findings propose a mechanism in which charge-dependent activation through ligand exchange enhances fibril recognition and promotes disruptive binding modes, demonstrating the potential of charge-tunable diruthenium complexes as therapeutic modulators of A{beta}42 fibril stability.
Hsu, I.-S.; Chou, Y.-C.; Lee, Y.-T.; Wang, W.-H.; Tsai, M.-Y.
Show abstract
Intrinsic tryptophan fluorescence is widely used as a sensitive reporter of protein conformational dynamics, yet the molecular origin of its temperature-dependent modulation remains unclear. Here we investigate the conformational dynamics of Trp134 in bovine serum albumin (BSA) using molecular dynamics (MD) simulations, free-energy calculations based on umbrella sampling and WHAM, quantum mechanical (QM) calculations, and QM/MM approaches. MD simulations show that the global structure of BSA remains stable while temperature induces a gradual population shift from the Ia+ to the Ia- rotamer. The corresponding free-energy landscapes reveal that this shift arises from subtle changes in basin stability and transition barriers along the rotameric coordinate. In contrast, standalone QM calculations on isolated tryptophan predict different energetic trends, highlighting the sensitivity of rotamer stability to electronic-structure treatments and environmental effects. QM/MM calculations partially reconcile these differences by incorporating the protein environment. Together, these results suggest that temperature reshapes the rotamer free-energy landscape of Trp134, leading to population shifts that modulate intrinsic tryptophan fluorescence in proteins.
Zajac, J. W. P.; Muralikrishnan, P.; Zeng, X.; Heldt, C. L.; Perry, S. L.; Sarupria, S.
Show abstract
Excipient effects on protein stability are critical for biological formulations, yet their selection remains largely empirical. Here, we use molecular dynamics simulations to define unifying metrics of protein-excipient interactions at atomistic resolution. Enhanced sampling simulations of fast-folding miniproteins, including Trpzip, WAAAH-helix (an alanine-rich -helix), and Trp-Cage, were performed to capture folding transitions across diverse excipient conditions. We identified a general stabilization mechanism based on shape complementarity between protein networks and surrounding solvent networks. Stabilizing excipients were found to form solvent structures that preferentially complement each protein, as well as residues central to known folding pathways. This framework enables a unifying approach to mechanism-based excipient selection across diverse protein and solvent chemistries. More broadly, by treating protein and solvent as dynamically coupled partners, it provides a transferable strategy for understanding solvent-mediated effects in complex molecular systems.
Weng, S. L.; Rekhi, S.; Kim, Y. C.; Palmer, J.; Mittal, J.
Show abstract
Biomolecular condensates exhibit spontaneous electrochemical microenvironments characterized by asymmetric ion distributions and pH gradients that emerge from protein-sequence-dependent charge regulation. Despite their biological importance, mechanistic understanding of these microenvironments has been constrained by the absence of computationally tractable frameworks capable of treating proton exchange, counterion partitioning, and buffer equilibria on consistent thermodynamic footing. Here, we introduce the buffered Charge-Regulation Monte Carlo (b-CR-MC) framework, which couples grand-canonical exchange of ions and buffer species with explicit charge regulation of titratable residues. By extending the CR-MC ion-merging strategy to multicomponent reservoirs and employing the Restricted Primitive Model, b-CR-MC achieves computational efficiency while maintaining thermodynamic rigor, with quantitative agreement to the more expensive generalized G-RxMC approach. Applied to full-length FUS (net positive) and PGL-3 (net negative) under physiological conditions, the framework reveals sequence-dependent pH gradients: the dense phase of FUS exhibits an alkaline shift, while PGL-3 exhibits an acidic shift, in both cases driving the condensate interior toward the protein's isoelectric point. Slab-geometry simulations further resolve the Donnan potential and continuous ion profiles across the condensate interface, confirming the direction and magnitude of these electrochemical shifts. Additionally, we identify spatially resolved buffer depletion within dense phases, establishing that dynamic charge regulation is a primary determinant rather than a secondary correction to condensate electrochemistry. By establishing a sequence-resolved, thermodynamically consistent computational platform, b-CR-MC enables quantitative prediction of how mutations and post-translational modifications reprogram condensate microenvironments across biological and pathophysiological contexts.
Pereira, A. F.; Araujo, J. O.; Tarraga, W.; Martinez, L.
Show abstract
Understanding the role of the protein backbone and side chains on cosolvent-induced stabilization is essential for a molecular picture of osmolyte action. The dominant view has been that protecting osmolytes stabilize proteins primarily through unfavorable interactions with the peptide backbone - the osmophobic effect - with side chains playing a minor or opposing role. By revisiting the decomposition of amino acid transfer free energies with proper account of the mutual shielding between backbone and side-chain groups, we derive a transfer model that is consistent with experimental denaturation m-values for urea and protecting osmolytes simultaneously - a feat neither the established nor the previously proposed universal-backbone models could achieve alone. A mechanism-dependent interpretation of backbone accessibility is proposed: geometric for excluded cosolvents, complete for binders where cosolvent-backbone interactions are specific. The model reveals that for all strong protecting osmolytes, including TMAO, sarcosine, sucrose, trehalose, and sorbitol, both backbone and side chains contribute favorably to protein stabilization, with side-chain contributions comparable to or exceeding those of the backbone. For urea, the model recovers the known balanced backbone and side-chain contributions to denaturation when the directional nature of urea-backbone hydrogen bonding is accounted for, which makes the backbone accessible to urea regardless of side-chain shielding. Weaker protectants such as proline, betaine, and glycerol are distinguished by competing backbone and side-chain effects that partially cancel. These results extend the osmophobic effect to protein side chains and establish a three-tier classification of osmolyte action: cooperative backbone and side-chain stabilization, cooperative destabilization, and competing contributions. The greater sensitivity of the model predictions to side-chain composition provides avenues for experimental validation of the underlying physical assumptions and for protein engineering.
Wieners, L.; Garcia, M. E.
Show abstract
The light absorption of the protein cryptochrome and its chromophore FAD is important for the regulation of circadian rhythms and in some species for sensing magnetic fields. To compute the absorption spectrum of chromophore, typically only a small region is treated quantum-mechanically due the high computational cost of spectroscopic calculations. We present a formalism that allows a quantum-mechanical treatment of not only the chromophore but also the neighbouring amino acids which differ from species to species. This is achieved by using the real-time time-dependent Hartree-Fock method. This method allows extending the quantum domain from typically only a few dozen atoms up to around 1,200 atoms for the largest calculations. The presented framework allows the treatment of neighbouring tryptophan residues or the cofactor molecule MTHF in the same calculation and allows to extract information of which regions absorb light depending on wavelength. The presented results also show that the environment around the chromophore FAD amplifies the light absorption in cryptochrome.
Majumder, A.; Dutta, M.; Cherek, L.; Voth, G. A.
Show abstract
HIV-1 buds from infected cells as immature virion particles with a scattered envelope glycoprotein (Env) distribution on their envelope. It then undergoes maturation, during which the viral protease cleaves the Gag polyprotein at multiple sites, leading to structural reorganization of the viral particle and lateral redistribution of Env proteins, ultimately rendering the virion infectious. However, the underlying mechanism of maturation-induced Env reorganization remains elusive. In this study, we combine microsecond-long all-atom (AA), bottom-up coarse-grained (CG) molecular dynamics simulations, and diffusion model-based backmapping to investigate the structural organization and key interactions of Env in viral membranes. AA simulations of fully glycosylated Env embedded in HIV-1 mimetic asymmetric bilayers were first performed to characterize its conformational dynamics and Env-lipid interactions. We then developed a bottom-up CG model of glycosylated Env from that AA data and simulated the mature HIV-1 virion envelope containing multiple Env proteins. The CG simulations predict that Env proteins form clusters through interactions mediated by the cytoplasmic tail domain (CTD) and adopt diverse tilted conformations within these clusters. These CG simulations were then backmapped to AA resolution and further AA simulations were carried out to identify, in detail, the specific interacting residues in the Env clusters. Additionally, analysis of epitope accessibility shows that broadly neutralizing antibodies (bnAbs) targeting the V1/V2 and V3 loops may efficiently interact with Env clusters on the mature virion surface. Together, these results provide a molecular mechanism for Env oligomerization during viral maturation and offer new insights into the accessibility of bnAb epitopes on Env clusters.
Mandal, M.; Shpinov, Y.; Lahlou, A.; Pham, F.; El Hajji, L.; Coghill, I.; Laureau, E.; Plamont, M.-A.; Perez, F.; Le Saux, T.; Aujard, I.; Gautier, A.; JULLIEN, L.
Show abstract
Reversibly photoswitchable fluorophores are widely used in advanced bioimaging but their design remains demanding. Here, we introduce a new series spanning the whole visible range, which results from combining a large set of fluorogens with the FAST protein scaffold. We first demonstrate that these well-established labeling fluorescent protein tags turn into negative reversible photoswitchers upon decreasing the fluorogen concentration and increasing light intensity. We then show that using not anymore one but two fluorogens adds new responses to illumination. Thus, we obtain positive reversible photoswitchers, that increase their brightness under illumination. We also generate a palette of non-covalent reversibly photoconvertible fluorescent proteins changing their fluorescence color upon illumination, a reversible behavior that still remains absent in regular fluorescent proteins. This light-induced color change opens the possibility to discriminate six spectrally similar FAST variants in live cells upon demonstrating the superiority of using multiple spectral channels for exploiting the time dependence of the fluorescence response to illumination.
Sen, A.; Chakrabarti, J.; Mitra, R. K.
Show abstract
The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.
Tanida, T.; Gofur, M. R.; Nakajima, T.
Show abstract
Forster resonance energy transfer (FRET) is a physicochemical phenomenon involving non-radiative energy transfer between donor and acceptor fluorophores. While FRET efficiency primarily depends on the proximity between fluorophores, additional factors also substantially influence the efficiency in living cells. However, how non-distance factors modulate live-cell FRET efficiency remains poorly understood. Here, we report the significant role of N- and C-terminal topology in determining live-cell FRET efficiency, independent of fluorophore proximity, donor variants, and subcellular compartment. Using acceptor photobleaching and sensitized emission measurements in living cells, we found that FRET efficiencies of mCherry-EGFP or mCherry-EYFP (acceptor-donor) were significantly higher than those of EGFP-mCherry or EYFP-mCherry (donor-acceptor), respectively. These efficiencies were comparable between the nucleus and cytoplasm. An orientation index analysis showed that the acceptor-donor configuration is more favorable than the donor-acceptor configuration regardless of donor variants and subcellular localization. FRET efficiencies were also higher with EYFP than with EGFP as the donor. AlphaFold2-based structural modeling suggested similar proximity with structurally heterogeneous and loosely constrained geometry of donor and acceptor fluorophores. Collectively, these results demonstrate that topological arrangement, rather than simple distance considerations, plays a significant role in FRET efficiency in living cells, providing molecular implications for the design of intramolecular FRET-based biosensors.